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)Q ABSTRACT 

Q\ The recent discovery by Cowie et al. and Tytler et al. of metals in the Lya clouds shows that 

the intergalactic medium (IGM) at high redshift is contaminated by the products of stars, and 
^ suggests that ionizing photons from massive star formation may be a significant contributor to the 
C/3 UV background radiation at early epochs. We assess the validity of the stellar photoionization 
\Q hypothesis. Based on recent computations of metal yields and O-star Lyman continuum (LyC) 
fluxes, we find that 0.2% of the rest-mass energy of the metals produced is radiated as LyC. By 
modeling the transfer of ionizing radiation through the IGM and the rate of chemical enrichment, we 
demonstrate that the background intensity of photons at 1 Ryd that accompanies the production of 
metals in the Lya forest clouds may be significant, approaching 0.5 X 10 -21 ergs cm -2 s _1 Hz -1 sr _1 
at z ~ 3 if the LyC escape fraction is > 0.25. Together with quasars, massive stars could then, 
If} in principle, provide the hydrogen and helium LyC photons required to ionize the universe at 
^ high redshifts. We propose that observations of the He II Gunn-Peterson effect and of the metal 
ionization states of the Lya forest and Lyman-limit absorbers should show the signature of a stellar 
spectrum. We also note that the stellar photoionization model fails if a large fraction of the UV 
radiation emitted from stars cannot escape into the IGM, as suggested by the recent Hopkins 
£^ Ultraviolet Telescope observations by Leitherer et al. of low-redshift starburst galaxies, or if most 
^ of the metals observed at z ~ 3 were produced at much earlier epochs. 
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1. INTRODUCTION 

It is widely believed that the background ultraviolet flux with energy greater than 1 Ryd 
arises from quasars (Sargent et al. 1980; Bergeron & Stasinska 1986; Madau 1992), possibly 
augmented by other sources such as hot, massive stars in star-bursting young galaxies (Bechtold et 
al. 1987; Steidel & Sargent 1989) or more exotic mechanisms such as decaying neutrinos (Sciama 
1993). This radiation is responsible for maintaining the ionization state of the diffuse intergalactic 
medium, the Lya forest clouds, and the metal-rich absorption systems. Using a model-independent 
argument from metal production, Songaila, Cowie, & Lilly (1990) suggested that the ionization 
state of intergalactic material could be attributed to star-forming galaxies if at least a few percent 
of galaxy formation occurred at redshifts z = 3 — 5. This possibility was studied in more detail by 
Miralda-Escude & Ostriker (1990). 

To date, a significant population of primeval, star-forming galaxies at high- z has not been 
detected. However, recent measurements (Tytler et al. 1995; Cowie et al. 1995; Fan 1994) show that 
a large fraction of the numerous, Lya forest clouds are associated with measurable C IV absorption 
lines. Because the Lya clouds have large filling factors and may contain a significant fraction of the 
baryons in the universe, the new data are the first evidence for a widespread distribution of metals 
in the IGM. If QSO absorption systems trace the bulk of star formation occurring in galaxies at 
high redshifts, and if the initial mass function (IMF) extends to massive O-type stars, the stellar 
contribution to the local ionizing background could be significant. 

In this paper, we examine the contribution of massive stars to the ionizing background, using 
a new calibration of the LyC emission from O stars, accurate radiative transfer of this radiation 
through the IGM, and a parameterization of the metal formation rate in QSO absorption-line 
systems. In § 2 we summarize the baryon content and metallicity of gas in Lya forest clouds, 
Lyman-limit systems, damped Lya systems, and star-forming galaxies. In § 3 we estimate the 
luminosity of ionizing photons that accompany the metal production and compute the spectrum 
of ionizing radiation. In § 4 we discuss the implications for the ionizing background and chemical 
evolution of the universe. 

2. HEAVY ELEMENT ABUNDANCES AT HIGH REDSHIFT 

2.1. Lyman-a Forest Clouds 

The term "Lya forest" is used to denote the numerous population of H I absorption systems 
observed in the spectrum of quasars (Weymann 1993; Bechtold 1994). With H I column densities 
ranging from 10 17 cm -2 down to 10 13 cm -2 and probably even lower, these systems are usually 
interpreted as optically thin clouds associated with the era of baryonic infall and galaxy formation. 
The clouds are observed to evolve rapidly with redshift between 1.7 < z < 4 (e.g., Murdoch et 
al. 1986) and typically do not show metal lines in ~ 1 A resolution spectra (Sargent et al. 1980; 
Tytler & Fan 1994; but see Lu 1991). However, recent spectra at high S/N and high resolution 
obtained with the HIRES spectrograph on the Keck telescope (Tytler et al. 1995; Cowie et al. 1995) 
have shown that 50%-60% of the Lya forest clouds with log iVni > 14.5 have undergone some 
chemical enrichment, as evidenced by weak, but measurable C IV absorption lines. The typical 
inferred metallicities, ZL ya , range from 0.003 to 0.01 of solar values, subject to uncertainties of 
photoionization models (e.g., Bergeron & Stasinska 1986; Steidel 1990; Madau 1991; Reimers & 
Vogel 1993) and possible contributions from local sources (Giroux, Sutherland, & Shull 1994; Viegas 
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& Gruenwald 1991; Gruenwald & Viegas 1993; Petitjean, Rauch, & Carswell 1994; Viegas & Friaca 
1995). Clearly, these metals were produced in stars that formed in denser gas; the metal-enriched 
gas was then expelled from the regions of star formation into the IGM. Although carbon need not 
be formed by massive stars, we assume that the IMF extends as a power law to O-type stars of 
mass 16 - 85 M©. 

The available data do not tightly constrain the mass density parameter, fiLya, °f the Lya 
clouds. Recent spectroscopic observations of close quasar pairs (Bechtold et al. 1994; Dinshaw et 
al. 1994) reveal that the forest clouds have characteristic sizes much larger than previously thought, 
with transverse radii > kpc, where h^o is the Hubble constant in units of 50 kms -1 Mpc -1 . 

Although the modeling is typically done for spheres, one can generalize by considering either spheres 
of radius R or clouds with large aspect ratios, a = R/£ = lOaio, of transverse size R to thickness £. 
Assuming photoionization equilibrium with the diffuse UV field and neglecting those unlikely cases 
in which the line of sight enters the cloud through its edge, we find that the total hydrogen gas 
density of a cloud of size R = R200 X 200 kpc and observed H I column density iVni = 10 14 iVi4 cm -2 
is 

1/2 

~ (5.6 x 10- 5 )R^ /2 Nll 2 Jig (a cos 6fl 2 cm" 3 , (1) 

where Thi is the hydrogen photoionization rate, 6 is the angle between the symmetry axis of the 
cloud and the line of sight, \ = n He/^H ~ 1/12, J912 = J-2i(z) x 10 -21 ergs cm -2 s _1 Hz -1 sr _1 is 
the metagalactic flux at the Lyman edge, «a = 4.1 X 10 -13 T 4 -0 ' 3 cm 3 s _1 is the (case-A) recombi- 
nation rate coefficient to all levels of hydrogen, and T = T4 X 10 4 K is the cloud temperature. For an 
ionizing radiation field with constant energy spectral index a, Thi ~ (1.2 X 10 -11 s _1 )J_2i/(a + 3). 
We adopt T4 = 4 and a = 0.5 for the background radiation. Note that setting (a cos 6) = 1 recovers 
the usual limit of spherical clouds. 

Press & Rybicki (1993) parameterized the observed frequency of absorbers per unit redshift 
by dN/dz ~ 4.2(1 + z) 2A6 for lines of rest equivalent width > 0.32 A, corresponding to columns 
Nm > 1-5 X 10 14 cm -2 (b = 30 km s _1 ). In an Einstein-deSitter cosmology (qo = 0.5), the 
volume filling factor of this subpopulation is / 4:.2(HoR/ a c cos#)(l + z) 4 ' 96 , assuming that 
all clouds have similar sizes. Denoting with ?^H,crit 

the closure hydrogen density, we find for 

^Lya = /^H/^H,crit at Z = 3, 

Q Lya * 0.02 h£ ( ZzH^V ' 2 . (2) 
\ a 10 cos 8 J 

We will show below that massive stars associated with QSO absorption line systems may generate 
an ionizing flux J-21 ~ 0.5 at these redshifts. 

In order to correctly estimate the contribution of the Lya forest to fl^ya over the entire 
range of column densities, a specific model for the dynamics of these clouds must be adopted. 
Values in excess of the cosmological baryon density of the universe derived from nucleosynthesis, 
O&/i 2 = 0.05 ±0.01 (Walker et al. 1991), would be obtained in the case of spherical clouds (Meiksin 
& Madau 1993; Press & Rybicki 1993). We conclude that the large Lya forest clouds must have 
large aspect ratios and contain a large fraction of the baryons in the universe. 
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The scarcer Lyman-limit systems (log A^hi > 17.2), which are optically thick to radiation with 
energy greater than 1 Ryd, are associated with the halo regions of bright galaxies, and represent 
essentially the same population of absorbers as the Mg II and C IV metal-line systems (Steidel 1993; 
Stengler-Larrea et al. 1995). Steidel (1990) used the column density measurements of H I and 
associated ions of heavy elements, in conjunction with photoionization models, to estimate the 
physical parameters of these systems. Although there are many uncertainties in the analysis, the 
typical metal abundances are Zlls ~ 0.003^0 at redshift z ~ 3, and their mass density parameter 
is 

LLS ~ O.OOS/^ 1 , (3) 

comparable to the mass density of visible matter in nearby galaxies. Recent observations obtained 
with the Keck telescope have confirmed that the metallicity of a sample of Lyman-limit systems 
with A^hi ~ 3 X 10 17 cm -2 is comparable to that of Lya forest clouds with A^hi ^ 3 X 10 14 cm -2 
(Cowie et al. 1995). 

2.3. Damped Lyman-a Absorption Systems 

The damped Lya systems (log A^hi > 20.3) may be the progenitors of present-day luminous 
galaxies (Wolfe 1990). While relatively rare, they account for most of the neutral hydrogen seen 
at high redshifts, with a mass density parameter evolving for 0.008 < z < 3.5 as 

n B (z) ~ (3.8 x 10" 4 )exp(0.83z)^ 1 , (4) 

for qo = 0.5 (Lanzetta, Wolfe, & Turnshek 1995). Thus, at z ~ 3, we have 0.005/i" 1 . At 

redshift z ~ 2, the typical metallicity, Zd, is approximately 1/10 of solar (Pettini et al. 1994). A 
"closed box" cosmic chemical evolution model which matches the observed evolution of Od(^) also 
implies that the metallicity of the damped Lya systems should decrease sharply at z > 3 (Lanzetta 
et al. 1995). 

2.4- Galaxies 

That star-forming galaxies might rival quasars as sources of ionizing radiation is easy to 
demonstrate. We assume, for purposes of estimation, that each population is dominated by objects 
near the break of their luminosity function. (This approximation is not rigorously true, since a 
galaxy luminosity function with a steep slope at L <C L* might result in dwarfs dominating the 
luminosity density.) Within factors of 3, the LyC photon luminosities of typical L* disk galaxies 
and quasars are given by S(L*) 10 53 s _1 and S(QSO) ~ 10 56 s _1 . For example, the total LyC 
production from the Milky Way was recently inferred by COBE observations of [N II ] 205 /xm 
emission to be 3.5 X 10 53 s _1 (Bennett et al. 1994), while the value quoted for QSOs is appropriate 
for a total luminosity ~ 10 46 ergs s _1 . If a fraction (f esc ) of the LyC from O stars escapes the 
galactic H I layer into the IGM, then the ratio of the LyC production rates of L* galaxies to that 
of QSOs at redshifts z = 2 — 4 is given by 



<p(QSO) 



S(L*) 
S(QSO) 



(/esc) ~ 10(/esc) , (5) 



where we have adopted comoving densities <f>(L*) = 10 2 Mpc 3 and <f>(QSO) = 10 6 Mpc 3 . 
Clearly, if (/ eS c) ~ 0.1, the two sources are comparable. From models of the escape of Galactic 
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LyC photons from OB associations through "H II chimneys" in the H I disk layer, Dove & Shull 
(1994) conclude that this escape fraction is at least 14%. The actual value could be somewhat 
higher if one accounts for the dynamical effects of superbubbles in puncturing the H I layer. We 
return to this issue in § 4. 

One problem with the above estimate is that we may not be justified in extrapolating the low- 
redshift density of bright, star-forming galaxies to z > 3. In fact, despite considerable observational 
effort, a significant population of primeval, young galaxies at high- z has to date evaded detection. 
While the bulk of the B < 24 galaxies are normal dwarfs at low redshifts (Colless et al. 1993), 
analyses of number counts, colors, and the lack of many Lyman-break candidates suggest that 
more than 90% of the faint, blue galaxies with 25 < B < 27 lie at z < 3 (Guhathakurta, Tyson, 
& Majewski 1990; Steidel & Hamilton 1993; Madau 1995). To reconcile these observations with 
the hypothesis that galaxies are major contributors to the ionizing background at high- z requires 
that we invoke dwarf or low surface brightness galaxies, undetectable to B = 27 at z > 3. Such 
galaxies are likely to have shallow gravitational potential wells, thus allowing easy ejection of heavy 
elements produced during their active star-forming phase. 

3. IONIZING PHOTON PRODUCTION 



3.1. The Ionizing Spectrum of Metal- Forming Galaxies 

If the metals detected in the Lya clouds arise from a population of stars with an IMF extending 
to massive stars (10 — 100 M©), this star formation will be accompanied by the production of hot 
gas, heavy elements, and ionizing stellar radiation. Sutherland & Shull (1995) used the evolutionary 
tracks of Schaller et al. (1992), the recently computed heavy element yields for massive stars and 
Type II supernovae (Woosley & Weaver 1995), and a recent recalculation of O-star LyC photon 
production rates (Vacca, Garmany, & Shull 1995) to derive a conversion efficiency of rest mass to 
ionizing continuum, 

^ y c ee J ° J " 2 « 0.002 . (6) 

Here, M m is the mass of heavy elements (Z > 6) produced in massive stars and L v is the spectral 
luminosity (ergs s _1 Hz -1 ) of ionizing photons with energies above hv = 1 Ryd. This conversion 
efficiency did not change by more than 10% when the IMF slope (Y) of stars on the upper main 
sequence [dN(> M)/dM = AM~ T for 8 M & < M < 85 M©] was varied from T = 1.0 to 2.0, where 
r = 1.35 corresponds to the usual Salpeter function (see also Songaila et al. 1990). The increased 
metal yields from high-mass stars are compensated for by a similar increase in the LyC production. 

Figure la shows the time-history of the ionizing spectrum from one such starburst. Except 
for the updated OB- and Wolf-Rayet stellar spectra, the resulting composite model is similar to 
the young-galaxy spectra of Bruzual & Chariot (1993). The stellar He II ionizing continuum above 
4 Rydbergs is dominated by Wolf-Rayet stars, for which we use the model spectra of Schmutz, 
Leitherer, & Gruenwald (1993). For the purposes of this paper, we shall assume that the starburst 
duration is sufficiently short compared to a Hubble time that the relevant quantity of interest is the 
time-integrated spectrum (Fig. lb). This should provide a reasonable estimate of the ensemble- 
average of a number of starburst populations. 



3.2. The Ultraviolet Background 
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The mean specific intensity J v of the radiation background at the observed frequency i/ Q , as 
seen by an observer at redshift z Q is: 

1 f°° dl (1 + z f 
J(u ,z ) = —J dz— ^ 1 + ^ 3 e(v,z)exp[-T e fi(v ,z ,z)] , (7) 

where e(v,z) is the proper volume emissivity (in ergs s _1 cm -3 Hz -1 ) at frequency v = v Q (l + 
z)/(l + z ), dl/dz is the line element in a Friedman cosmology, and exp[— r e ff (i/ Q , z Q} z)] is the 
average transmission of a clumpy universe. We have computed the effective photoelectric optical 
depth r e ff = — ln((e _T )) of a cloudy intergalactic medium following Madau (1995). For Poisson- 
distributed clouds we have (Paresce, McKee, & Bowyer 1980) 

f z f°° d 2 N 
T ef£ (v 0} z 0} z) = / dz' / dN m — — — — (1 - e T ) , (8) 
Jz Jo dN m dz' 

where d 2 N/ dNjudz' is the redshift and column density distribution of absorbers along the line 
of sight, and r is the H I continuum optical depth through an individual cloud of column density 
Nhi- When r <C 1, T e fi becomes equal to the mean optical depth. In the opposite limit, the 
obscuration is picket fence-type, and the effective optical depth becomes equal to the mean number 
of optically thick systems along the line of sight. The cosmic LyC opacity is dominated by absorbers 
with t ~ 1. Note that the He I contribution to the attenuation is negligible if QSOs contribute 
substantially to the ionizing background, while He II absorption 1 on the way must be included 
when (1 + z) > (1 + z )(X /228 A). 

The ionizing emissivity associated with metal formation can be written as 

t(y,z) = joo^ ^ VLyC O abs /9crit(^)c 2 Z(/ esc ) , (9) 

where the first term in parentheses is the normalized, time-integrated spectrum of our template 
star-forming galaxy (see Fig. 16), O a b s is the baryonic mass density parameter of metal-forming 
material in QSO absorption-line systems, and Z is the metal formation rate. To parameterize in 
simple form the many possible evolutionary models for the history of cosmic metal formation, we 
assume that 

Z = const = - — - — , ^ °^ . — , (10) 

[Mz ) -< H (Zon)] 

where £h is the Hubble time, and z on flags the cosmic epoch at which metal formation turns on. 
Equations (7), (9), and (10) yield a background flux proportional to the product a b s -^(^ o )(/esc)- 
As discussed in § 2, the gaseous baryons in the IGM add up to O a b s in the range 0.01h^ to 0.05h^ , 



1 An approximate (within 5%) formula for the effective H I continuum opacity at the ob- 
served wavelength A G < 912 A is r eff (A ,z ,z) ~ 0.25x 3 (x 0A6 - x o A6 ) +9.4^- 5 (x - 18 - x o 18 ) 
-Q.lxlix- 1 ^ 2 - x" 1 - 32 ) -0.023(x 1 - 68 -x 1 - 68 ), where x c = (A G /912 A)x Q , x Q = (1 + z Q ), and 
x = (1 + z) (Madau 1995). The first term on the right-hand side represents the contribution 
of Lya clouds, the others are due to Lyman-limit systems. 
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while the metallicities Z range between 0.01 and 0.1 times the solar value of 0.02. Although the 
product fi a bs.Z may be highest in the damped Lya systems, the LyC escape fraction from such 
systems may be lower. Note that, at high redshifts, the ionizing background will depend only 
weakly on the assumed metallicity evolution history. At early epochs, sufficient neutral hydrogen 
is contained in the IGM to significantly attenuate the LyC flux from all but the nearby sources, 
within Az 0.2. Because of this large optical depth, our estimate of the metagalactic flux at 
z Q 3 will depend only on the "local" metal formation rate Z. 

Figure 2 shows the 300 — 1000 A spectrum of the metagalactic radiation field from massive 
stars, computed at z Q = 3 for different values of z on = 3.2, 3.5, and 4.0 and filtered through a 
cloudy IGM. We assume a value of Q !l \ )S Z(z )(f esc ) equal to 1 X 10 -6 , corresponding for example 
to fi a bs = 0.02 (comparable to the baryon density in the Lya forest clouds), Z = 2 X 10 -4 (about 
one percent of solar, the inferred metallicity of the Lya forest and Lyman-limit absorbers), and 
(/esc) = 0.25. This product is quite uncertain; we return in § 4 to a discussion of values of the 
individual terms. In the three cases considered, the background intensity at the hydrogen Lyman 
edge, J912, is found to range from 0.13 to 0.36 X 10~ 21 hl ergs cm -2 s _1 Hz -1 sr _1 . These values are 
a substantial fraction of direct measurements of the metagalactic ionizing flux from the "proximity 
effect" (Bajtlik, Duncan, & Ostriker 1988; Bechtold 1994) and are comparable to the estimated 
contribution from observed QSOs (Meiksin & Madau 1993). 

4. DISCUSSION 

Recent observations of the cosmological mass density of metals associated with QSO absorption 
systems at high- z are expected to bear significantly on the history of star formation and the chemical 
and ionization evolution of the universe as a whole. We have shown that the LyC photons which 
accompany the production of metals in the Lya forest clouds could contribute significantly to 
the ionizing background at high redshifts. In this scenario, the metals are formed from massive 
stars, probably in dense gaseous regions of young galaxies, and are then expelled to large distances 
consistent with the sizes of Lya clouds. As a consequence, it is difficult to ascertain the exact 
location of the LyC production except in an average sense. Similarly, the photoionization rates 
derived here do not distinguish between "local" sources (from a single star-forming galaxy) and 
the integrated effects of many such galaxies. One would naturally expect a detailed model to 
include fluctuations in the spectrum due to varying local star-formation rates and "filtering" by 
surrounding H I gas (the parameters Z and (f eS c) in eq. [9]). 

The model calculations we just described are based on a number of parameters, some of 
which contain significant uncertainties. In order to provide an assessment of the validity of the 
stellar photoionization hypothesis, we believe that it is important to summarize the weak spots 
in the argument. As shown in equation (9), the LyC emissivity is proportional to the product 
QnhsZ(z )(f esc ) , multiplied by a spectral template and an efficiency, ??LyC ~ 0.002, with which 
LyC radiation is generated during the production of heavy elements. The spectral template is a 
fairly robust quantity, since it is based on the composite, time-integrated spectrum of massive stars 
produced in a burst of star formation. Likewise, the efficiency ??LyC is nearly constant, as long as 
the stellar population extends as a power-law to massive stars (50-80 M©). The metallicity of 
the Lya clouds was based on detection of absorption by C IV , and not oxygen. Since carbon is 
produced primarily by low-mass stars, there is always the possibility that star formation at high 
redshift may have occurred with an anomalous IMF lacking massive stars. We find this possibility 
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slim, but the detection of VI in the Lya forest clouds would settle the issue, since oxygen is 
unambiguously produced in massive stars (Woosley & Weaver 1995). 

Next, we address the possible uncertainties associated with our estimate of the starburst 
contribution to the ionizing background. We choose to relate these issues to uncertainties in the 
three key terms in equation (9): the baryon density, fi a bs, the metallicity, Z(z Q ) } and the LyC 
escape fraction, (f eS c)- The first term is on a fairly solid basis; the baryon density associated with 
all types of Lya absorbers at z ~ 3 is at least 0.01 of closure density and could approach the 
nucleosynthesis limit, (0.05 ± 0.01)/i^~ 2 . The Lya forest potentially contains the bulk of these 

baryons, although the uncertainties due to cloud ionization and geometry are significant. 

The second term, defining the cloud metallicity history, is somewhat more uncertain. Although 
the metallicities at a given epoch are as reliable as the ionization models used to derive them, the 
time history of metal production is less well defined. Current estimates for damped Lya absorbers 
(Pettini et al. 1994) suggest that the bulk of metals formed between redshifts z = 2 — 3; no such 
estimates exist for the Lya forest. For example, the metal formation rate near the observer redshift 
z Q could be negligible, and most of the element production may occur at earlier epochs. In this 
case, the background field at z Q would be quenched by the large continuum optical depth of the 
universe at those redshifts. More likely, galactic chemical history, like star formation history, is 
spread widely over time. In this case, our approximation of a constant metal formation rate may 
not be too bad. 

The third term, (/ esc ), parametrizing the LyC escape fraction, is probably the most uncertain 
of all. This factor depends on a variety of stellar, galactic, and interstellar quantities, with few 
observations to guide us. In a theoretical study of the structure of H II regions in vertically stratified 
galactic H I disk layers punctured by OB associations, Dove & Shull (1994) used distributions of 
H I column density (Dickey & Lockman 1990) appropriate for the solar vicinity, together with the 
ionizing photon luminosity function of OB associations (Kennicutt, Edgar, & Hodge 1989). They 
concluded that 14% of the LyC photons escaped through the H I layers into the halo in the 2.5 
kpc region around the solar circle, and they speculated on how this fraction might change in other 
locales. For example, the escape of LyC from a large star-forming complex in a galactic nucleus 
would likely be higher than that from a distribution of smaller OB associations in a galactic disk. 
Similarly, the gas distributions and OB formation sites in dwarf or irregular galaxies probably differ 
from those in spiral disks. Finally, there is almost no information characterizing the structure and 
types of star-forming galaxies at high redshift. 

A few observational tests might provide some help in constraining (f eS c)- Using spectral 
observations with the Hopkins Ultraviolet Telescope, Leitherer et al. (1995) have recently set limits 
of a few percent on the LyC escaping fraction in low-redshift starburst galaxies. Thus, at first sight, 
the stellar photoionization hypothesis may be untenable, at least at low redshift. However, even 
though their target galaxies were chosen to be bright in the near-ultraviolet, there is no guarantee 
that the LyC is not obscured by gas along a different orientation (cone effects). In addition, these 
constraints need to be obtained for dwarf galaxies as well. Searching for Ha or Lya emission might 
also constrain the escaping flux, although the quantitative sensitivity of this test is likely to be 
diminished by internal dust absorption (Binette et al. 1993). The interpretation of Lya emission 
and LyC radiation escaping from galaxies is intrinsically difficult, since they rely on modeling of 
optical and near-UV radiation to determine the total ionizing radiation produced. We conclude, 
therefore, that (f eS c) is an unconstrained parameter of our model and of previous models (Songaila 



- 9 - 



et al. 1990; Miralda-Escude & Ostriker 1990) of stellar contributions to the intergalactic radiation 
field. 

Despite these caveats, we believe that our estimates provide some rationale for the hypothesis 
that hot stars could rival quasars as a source of photoionization of the IGM. If this is true, the 
signature of a stellar spectrum should appear in several places, notably the He II Gunn-Peterson 
effect, whose detected strength (Jakobsen et al. 1994) requires a relatively soft radiation field 
at 4 Ryd compared to 1 Ryd (Madau & Meiksin 1994; Giroux, Fardal, & Shull 1995). On the 
other hand, if N(He II )/N(H I ) in the Lya clouds is measured to be less than 100, then stars are 
unlikely to be major contributors to the ionizing background; after filtering by the IGM the stellar 
spectra are too steep. Likewise, the distribution of metal ionization states of the Lya forest and 
Lyman-limit absorbers should reflect an ionizing stellar contribution, from either local sources or 
the integrated background. 

Therefore, we propose three crucial tests of the stellar photoionization scenario: (1) confirming 
the presence of massive stars in the high- z stellar IMF by detecting O VI AA1032, 1038 absorption 
in the Lya clouds; and (2) constraining the LyC escape fraction, (/ eS c) 5 by far-UV continuum 
and Lya emission observations of low-redshift starburst galaxies and by indirect estimates for star 
forming regions in the Galaxy; and (3) measuring the ratio N(He II )/N(H I ) in the lines of sight 
to high- z quasars. 

We thank Ralph Sutherland for his help in deriving the "starburst galaxy" ionizing spectrum, 
and we acknowledge helpful conversations on the subject of this paper with Mark Giroux, John 
Stocke, and David Tytler. Support for this work was provided by the Colorado Astrophysical The- 
ory Program (NASA grant NAGW-766). JMS thanks the Visitor Program at the Space Telescope 
Science Institute where this work was completed. 
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FIGURE CAPTIONS 



Figure 1: Modeled input spectrum (Sutherland & Shull 1995) from a gaussian starburst of width 
2cr = 4 Myr containing 5000 Mq of stars between 8 to 85 Mq with IMF slope T = 1.6 (see 
§ 3.1). Spectrum shows four bands: soft UV (5.0 - 13.6 eV) and ionizing continua of H I (13.6 
eV), He I (24.58 eV), and He II (54.4 eV). (a) Time history of ionizing spectral luminosity at times 
2, 4, 6, 12, and 16 Myr (peak is at 6 Myr). (b) Time-integrated spectrum (E v in eq. [9]) of the 
starburst normalized to the total radiated ionizing energy E v dv , integrated over the 42.5 Myr 
lifetime of an 8 Mq star. 



Figure 2: Observed spectrum at redshift z Q = 3 of a UV background dominated by massive stars, 
modified by the absorption of intervening Lya clouds and Lyman-limit systems. We assume a 
constant metal-formation rate from z on to z Q and adopt Q !l \ )S Z(z )(f esc ) = 1 X 10 -6 (see § 3.2). 
The metagalactic flux is plotted in units of (hf ) ergs cm -2 s _1 Hz -1 sr _1 . From top to bottom, the 
three curves correspond to different epochs, z on = 3.2, 3.6, and 4.0, for the onset of heavy element 
formation. The spectrum for z on = 3.2 shows bumpy structure because the stellar input spectrum 
undergoes less redshift smearing. 
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